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ABSTRACT 

We present a detailed theoretical analysis of the gravitational- wave (GW) signal of the post-bounce 
evolution of core-collapse supernovae (SNe), employing for the first time relativistic, two-dimensional 
(2D) explosion models with multi-group, three-flavor neutrino transport based on the ray-by-ray-plus 
approximation. The waveforms reflect the accelerated mass motions associated with the characteris- 
tic evolutionary stages that were also identifled in previous works; A quasi-periodic modulation by 
prompt postshock convection is followed by a phase of relative quiescence before growing amplitudes 
signal violent hydrodynamical activity due to convection and the standing accretion shock instability 
during the accretion period of the stalled shock. Finally, a high-frequency, low-amplitude variation 
from proto-neutron star (PNS) convection below the neutrinosphere appears superimposed on the 
low-frequency trend associated with the aspherical expansion of the SN shock after the onset of the 
explosion. Relativistic effects in combination with detailed neutrino transport are shown to be essen- 
tial for quantitative predictions of the GW frequency evolution and energy spectrum, because they 
determine the structure of the PNS surface layer and its characteristic g-mode frequency. Burst-like 
high-frequency activity phases, correlated with sudden luminosity increase and spectral hardening of 
electron (anti-)neutrino emission for some 10 ms, are discovered as new features after the onset of the 
explosion. They correspond to intermittent episodes of anisotropic accretion by the PNS in the case 
of fallback SNe. We find stronger signals for more massive progenitors with large accretion rates. The 
typical frequencies are higher for massive PNSs, though the time-integrated spectrum also strongly 
depends on the model dynamics. 

Subject headings: supernovae: general — neutrinos — radiative transfer — hydrodynamics — gravitation- 
gravitational waves 
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1. INTRODUCTION 

Core-collapse supernovae have been considered as a 
source of gravitational waves (GWs) since the 1960s 
(IWeberi [T966'l. While they may have been ousted by 
compact binary mergers as the most promising source for 
the first direct detection of GWs by now, the prospec- 
tive GW signal from a nearby supernova may still pro- 
vide enormously important clues about the dynamics 
in the supernova core, shed light on the nature of the 
explosion mechanism, and could possibly provide con- 
straints for the equation of state (EoS) of neutron star 
matter. If supernovae are to become a fruitful subject 
of the future field of observational GW astronomy, this 
will require reliable waveform predictions, and these can 
only be based on simulations that accurately capture the 
evolution from the collapse through several hundreds of 
milliseconds of accretion by the stalled shock front into 
the explosion phase - a formidable challenge considering 
the host of different factors (neutrino transport, multi- 
dimensional hydrodynamical instabilities, general rela- 
tivity, EoS physics, etc.) that infiuence the dynamics. 

There is a number of possible scenarios for GW emis- 
sion from supernovae. Rotating progenitors already pro- 
duce a signal during the phases of collapse and bounce 
due to the breaking of spherical symmetry, a scenario 
which has long received interest from the numeri cal rel- 
ativity community in particular (see lOttl 120091 for an 
extensive review). State-of-the art predictions of the 
GW signal come from 2D and 3D general relativistic 



simulations (lOtt et al.ll2q07al: iDimmelmeier et al.ll2007al 



[2008t lAbdikamalov et al.ll2010l: lOtt et al.ll2012[ ) emplov- 
ing a parametrized "deleptonization s cheme" for the 
core-collapse phase (jLiebendorfeii 120051 ) . and, most re- 
cently, a neutrino leakage treatment for the post-bounce 
phase (Ott et al. 2013). The accuracy of this current ap- 
proach is yet to be tested against self-consistent models 
including neutrino transport, but it is conceivable that 
the dynamics of rotational collapse can be captured rea- 
sonably well even with a simplified neutrino treatment. 
However, the fact that most sup ernova progenitor s are 
believed to rotate rather slowly (jHeger et al.ll2005[) im- 
plies that the signal from rotational core bounce may be 
weak and difficult to detect. 

The situation is different for the GWs produced 
by the different hydrodynamical instabilities that de- 
velop during the post-bounce phase also in non- 
rotating progenitors, such as convection in the 
proto -neutron star (as first pointed out by 'Eostein 



1979 ). in the n eutrino-hea t ed ho t- bubble region (iBethe 
Tggg Hcrant eFaLl IT991 [1991 lBurrows et~aII 119951: 



Janka fc M(illedll99a IMiiller fc Jankalll997D. the stand- 



ing accretion-shock instability ("SASI", Blondin et al.l 
2003; Blondin & Mczzacaooa 2006; Foglizzo et al. "20061 
Ohnishi et al. 2006: Fogl izzo et al. 2007; Scheck et al| 
2008; Iwakami et al. 20081 l2009HFernandez fc Thompson! 
2009i; ,Fernandez .2010'). 



Regardless of the actual nature of the explosion mech- 
anism, the evolution during the post-bounce phase de- 
pends crucially on the effects of neutrino heating and 
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cooling and thus requires an elaborate treatment of 
the microphysics and in particular of the neutrino 
transport. Because of this constraint, studies ad- 
dressing the GW signal from the first several hun- 
dreds of milliseconds of the post-bounce evolution have 
been limited either to the Newtonian approximation 
or the "pseudo-New tonian effective potential" approach 
(|Marek et all |2006[) until now. Since the 1990s, sev- 
eral authors have investigated the problem in 2D and 
3D, mostly relying on parametrized approximations for 
the neutrino heating and cooling (Miillcr & Janka 1997,; 
Kotake et all l2007t iMurphVet al. 2009; Kotake_eLaD 



2009D, the IDS A method (iScheidegger et al. ,2010l). or 
on gray neutrino tr ansport schemes ( Frver et al.l |200'S 
iMiiller et al.l 120 12H ). In addition, gravitational wave- 
forms from state-of-the art multi-group neutrino hydro- 
dynamics simulations of core-collapse supernova e have 



becorne available dur i ng the recent years (Miiller et al 



200l lOtt et al.|[200l [Marek et al.l 120091: lYakunin et al 



2010| ). These different (pseudo-)Newtonian studies have 



by now established the qualitative features of the GW 
signal from the post-bounce phase: During the first 
several tens of milliseconds, prompt post-shock convec- 
tion gives rise t o a quasi-periodic signal in the range 
around 100 Hz (iMarek et al.| [200l iMurphv et all [20091: 
lYakunin et al.l [20101 ) . which is followed by a period of 
reduced GW activity until hot-bubble convection and 
the SASI become vigorous and lead to the emission of 
a stochastic sign al with typical frequencies of several 
hundreds of Hz (IMiiller et al.l [200l IMarek et al.l [200l 
IMurphv et al.ll2009[ ). After the onset of the explosion, 
asymmetric shock expansion may produc e a "tail", i.e. 
a growing offset, in the GW amplitude (IMurphv et al.l 
[2001 lYakunin et all [20T0I : IMiiller et all l2012d ). and~a 
high-frequency signal from proto-neutron star convec- 
tion starts to appear. Tail-like sign als may also result 
from anisotropic n e utrino emission (IMiiller et al.l[200l 
Marek et all 120091: lYakunin et all l2010t IMiiller et al l 



2012c[ ) powered by accretion downflows onto the proto- 



neutron star. 

However, the accuracy of the Newtonian approxima- 
tion and even of the "effective potential" approach, 
which have formed the basis of GW predictions for the 
post -bounce phase so far, is limited. Bo t h ID stud- 
ies (iBaron et all 119891: iBruenn et all 120011 : iLentz et al.l 
[20TI" as w ell as the recent 2D models of our own group 
(jMiiller et al. 2012b) have demonstrated that GR effects 
have a non-negligible impact on the neutrino emission, 
the shock position, and the heating conditions in the 
supernova core. It is conceivable that GR affects the 
GW signal to a similar degree. To address this question, 
we present gravitational waveforms from axisymmetric 
(2D) general re lativistic simulations (using t he xCFC ap- 
proximation of iCordero-Carrion et al.l[2009f ) of the post- 
bounce phase including multi-group neutrino transport 
for the first time. We compare our results with pre- 
dictions from Newtonian and pseudo-Newtonian mod- 
els. We discuss the relativistic GW signals of six dif- 
ferent progenitors with zero-age main sequence masses 
ranging from 8.1Mq to 27Mq (among them five ex- 
plosio n models), part of w hich have alread y been stud- 
ied in IMiiller et al.l ()2012b[ ) (paper II) and iMiiller et al] 
(|2012air " For a 15Mq star, two simulations using the 
Newtonian approximation and the effective potential ap- 



proach provide the basis for diagnosing GR effects and for 
working out the reason of systematic differences where 
possible. We also consider a ISMf^ model with sli ghtly 
simplified neutrino rates from IMiiller et all (|2012bl ). 

Our paper is structured as follows: In Section [21 we 
briefiy summarize the numerical methods a nd the model 
setup, which is laid out in greater detail in IMiiller et al.l 
(|2010l I2012bl paper I and II). In Section [3l we describe 
the methods used for extracting and analyzing the grav- 
itational wave signals. The wave signals produced dur- 
ing the different phases of the post-bounce evolution - 
the early phase of shock propagation and prompt post- 
shock convection, the steady-state accretion phase, and 
the explosion phase - are discussed in Section [4| on the 
basis of two exemplary explosion models. The depen- 
dence of the gravitational wave emission on the progen- 
itor model is discussed separately in Section [5l Finally, 
we summarize the most salient features of our relativis- 
tic waveforms in Section [Hi Important technical issues 
relied upon in our paper, such as the derivation of a 
modified GW quadrupole formula, analytic expressions 
for the GW signal produced by an aspherical shock front, 
and the Brunt- Vaisala (buoyancy) frequency in GR are 
treated in Appendices A, B and C. 

2. NUMERICAL METHODS AND MODEL SETUP 

We analyze the GW emission for several 2D 
simulations performed with the general relativis- 
tic neutrino h ydrodynamics code Vertex-CoCoNuT 
(jMiiller et al.l [2010.) or its (ps eudo)-Newtonian coun- 
terpart Vertex- Prometheus (jRampp fc Jankal 120021 : 
iBuras et al.l 120061 ). The reader should refer to these pa- 
pers for an in-depth description of the numerical methods 
and the input physics (including the choice of progenitor 
models, neutrino interaction rates, and the equation of 
state). At this point, we confine ourselves to a summary 
of those technical aspects that are directly relevant for 
the present paper with its focus on GW signals. 

The hydrodynamics modules CoCoNuT and 
Prometheus used in conjunction with the neu- 
trino transport solver Vertex both rely on similar 
higher-order finite-volume techniques, but differ in their 
treatment of general relativity. Prometheus solves 
the equation of Newtonian hydrodynamics, but can 
optionally account for some effects of strong-field gravity 
by means of a modified gravitational potential ( "effective 
potential"; iRampp fc Janka 2002; Mare k et al. 2006^ ). 
By contrast, the general relativistic equa tions of hydro- 
dynam i cs are solved in CoCoNuT (Dimmelmeier et al.l 
l2002al 120051 ). and th e extended conformal flatn ess 
approximation (xCFC, ICordero-Carrion et al.l 120091 ) is 
used for the space-time metric. Comparisons of xCFC 
with the full ADM formalism in the context of rotational 
core-c oll apse have shown excellent ag reement (jOtt et al.l 
I2007al[bl : iDimmelmeier et al.l l2007bl) and suggest that 
it is fully adequate for capturing the dynamics in the 
supernova core. However, xCFC is a waveless approxi- 
mation so that we need to extract the GW signal in a 
post-processing step with the help of a suitable version 
of the quadrupole formula (see Section [3]) . 

The neutrino transport module Vertex solves the 
energy-dependent neutrino moment equations for all neu- 
trino flavors using a variable Eddington factor technique 
(fRampp fc Janka. 2002. ) and relies on the so-called "ray- 
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Figure 1. Entropy along the north and south polar axis as a function of time for the relativistic simulations G8.1, G9.6, Gil, G15, 
G25, and G27 (top left to bottom right). The shock trajectory is visible as a discontinuity between the darker violet and blue tones in the 
pre-shock region and lighter colors (red, green, yellow) in the post-shock region. 



by-ray-plus" ap proximation for m ulti-dimensional neu- 
trino transport (jBuras et al.ll2006() . The ray- by-ray-plus 
approach allows us to predict angular variations in the 
neutrino radiation field (and hence the low-frequency 
GW signal generated by the neutrinos) at least in rough 
qualitative agreement with full rn ulti-angle transport 
(jptt ct al. 2008; Brand t et al.ll2011[ ). 

In total, we simulate the evolution of six different pro- 
genitor models with zero-age main sequence masses rang- 
ing from S.IM© to 27Mq (simulations G8.1-G27). These 
include a metal-poor S .lMp) (10~'^ solar rn etallicity) pro- 
genitor (model uS.l in lMiiller et al.|[2012aD . and a metal- 
free 9.6M0 star (z9.6, Alexander Heger, private com- 
munication), while the other progeni tors (models sll.2, 
s25.0 and s27.0 of Woos ley et al.|[2002l and model sl5s7b2 
of iWooslev fc Weaver., 19951) have solar metallicity. Three 
of the progenitors (z9.6, s25. 0, s27.0) were simula t ed us - 
ing the equation of state of 'Latti mer fc Swestvl (|1991[ ) 
with a bulk incompressibility modulus of nuclear mat- 
ter oi K ^ 220 MeV (LS220), whereas K = 180 MeV 



(LS180) was applied in all other cases. Because of 
very similar proto-neutron star radii for the neutron star 
masses encountered in this study, LS220 and LS180 yield 
very similar results, and the use of LS180 for some of the 
progenitors is justified de spite its marginal inco nsistency 
with the I.97M0 pulsar (i Demorest et al.ll2010D . For de- 
tails, we refer the reader to Paper II. 

The general features of the gravitational wave signal 
from the post-bounce accretion and explosion phases 
are discussed on the basis of the II.2M0 progenitor of 
'Wooslcv et al.' (2002^ and the ISM© model sl5s7b2 of 
Woosley & Weaver (1995]) as prototypes for "early" and 
"late" explosions. For the 15Mq progenitor, we also con- 
sider three models with a different treatment of gravity 
(G15: GR; M15: effective potential; N15: purely New- 
tonian), as well as a GR model (S15) with simplified 
neutrino rates in order to explore the impact of these 
factors on the GW signal. The simplifications in model 
S15 include t he use of the FFN rates (Fuller ct al. 198^ 
lBruennlll985D for electron captures on nuclei, a simpler 
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treatment of neutrino-nucleon reactions, and the omis- 
sion of neutrino-neutrino pair conversion (see paper II 
for details). 

Except for the 25Mq case, explosions have been ob- 
tained for all progenitors with GR and with the full set 
of neutrino rates. Figure [T] provides a compact overview 
over the relativistic (G-)series of simulations: Models 
08. 1 and G9.6 explode rather early and exhibit convec- 
tive activity only on a moderate level after the onset of 
the explosion. The 11.2Mq model Gil shows a much 
slower expansion of the shock and several violent shock 
oscillations before the explosion takes off. Model G15 de- 
velops a very asymmetric explosion as late as ~ 450 ms. 
The more massive 25Mq and 27Mq models G25 and G27 
differ from the other models by a more clearly discernible 
SASI activity, visible as strong periodic sloshing motions 
of the shock in Figure [U which lead to an explosion in 
the case of G27. We note that no explosion develops 
in the simulations without GR and/or the full neutrino 
rates (M15, N15, S15). 

A summary of all nine models considered in this paper 
is given in Table [TJ For a detailed discussion of models 
G11.2 and G15, sec Mii ller ct al. (2 12b), a nd for details 
on G8.1 and G27, see l'MiiUer et all (|20T2l . 

3. GRAVITATIONAL WAVE EXTRACTION 

The xCFC approximation used in Vertex- 
CoCoNuT does not allow for a direct calculation 
of gravitational waves as the corresponding degrees 
of freedom in the metric are missing. We therefore 
need to extract gravitational waves in a post-processing 
step with the help of some varia nt of the Einstein 
quadrupole formula ()Einsteinlll918() . Modified versions 
of the Newtonian quadrupole formula (exploiting am- 
biguities concerning the identification of Newtonian 
and relativistic hydrodynamical variables) have been 
found to be reasonably accurate eve n in the strong-field 
regime JShibata & Sekiguchil [20031 : INagar et al\\2(M : 
iCordero^ arrion et al. 2012)^ For the gauge used in 
Vertex-CoCoNuT and the typical conditions in a 
supernova core, it is possible to derive a modified version 
of the time - integrated Newtonian quadrupole formula 
(iFhml [TqsI iFinn fc Evailsl \Tml IBlanchet et aD [iQQOh 
directly from the field equations (see Appendix |A|) . 
Assuming axisymmetry, we obtain the quadrupole 
amplitude A^q in non-geometrized units as 
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Here, (/) is the conformal factor for the three-metric in the 
CFG spacetime, and Si denotes the components of the 
relativistic three-momentum density in the 3-1-1 formal- 
ism, which is given in terms of the rest-mass density p, 
the specific internal energy e, the pressure P, the Lorentz 
factor W ^ and the covariant three- velocity components Vi 
as 



S,=p{l + e + P/p)W\, 



(2) 



Si^u denotes the momentum source term for Si due to 
neutrino interactions (which must be subtracted from 
dSi jdt as explained in Appendix |X| . In practice, these 
neutrino source terms do not yield a significant contri- 
bution to the integral in Equation ([1]). 

Afp determines the dimensionless strain measured by 
an observer at a distance i? and at an inclinatio n angle 
9 wit h respect to the equatorial plane (see, e.g.. iMiilleii 

[Till), 
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In the following, we will always assume the most opti- 
mistic case of an observer located in the equatorial plane, 
i.e. sin^ 9—1. In addition to the gravitational wave sig- 
nal from the matter, we compute the gravitational wave 
signal due to anisotropic neutrino emission using the Ep- 
stein formula (Epstein 1978; M idler &: Janka.l997i) for 
the gravitational wave strain h^, 
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L,{t')a,{t')dt' 



(4) 



Here is the total angle-integrated neutrino energy 
flux, and the anisotropy parameter can be obtained 

8-8 \ r dL 

= — / 7rsin6'(2|cos0| - 1) — fdr2 (5) 
i-iv J dsZ 

in axisymmetry (jKotake et al.l 120071 ) . hy can be con- 
verted into an amplitude A^q ^ by inverting Equation 

The energy £^gw radiated in gravitational wa ves can 
be co mputed from A^q ^ as follows (see, e.g., iMulleii 

[Tm, 
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dt 



We also calculate the spectral energy distribution dE/df 
of the gravitational waves, 
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(7) 



Previous studies ([Marek et ahl l2009t iMurphv et"all 
120091 iMiiller et all l2012cf) have relied on Fourier trans- 
form methods (such as the short-time Fourier transform) 
for analyzing temporal variations in the frequency struc- 
ture of the signal. While the Fourier transform of the 
signal can be related directly to the power radiated in 
gravitational waves, the time-variable frequency struc- 
ture of the signal can be captured more sharply with 
wavelet transforms. 

We therefore compute the wavelet transform x(i, /) 
(expressed as a function of time and frequency) of the 
gravitational wave signal using the Morlet wavelet (see, 
e.g.. E quations 1 and 2 and Table 1 of lTorrence fc Compol 
Il998f) with wavenumber fc = 20 and define a normalized 
power spectrum x(i, /): 
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Table 1 

Model setup 







neutrino 


treatment of 


final post-bounce 


angular 


explosion 


time of 




model 


progenitor 


opacities 


relativity 


time reached 


resolution 


obtained 


explosion'' 


EoS 


G8.1 


u8.1 


full set 


GR hydro - 


h xCFC 


325 ms 


1.4° 


yes 


175 ms 


LS180 


G9.6 


z9.6 


full set 


GR hydro - 


h xCFC 


735 ms 


1.4° 


yes 


125 ms 


LS220 


Gil 


sll.2 


full set 


GR hydro - 


h xCFC 


750 ms 


2.8° 


yes 


213 ms 


LS180 


G15 


sl5s7b2 


full set 


GR hydro - 


h xCFC 


775 ms 


2.8° 


yes 


569 ms 


LS180 


S15 


sl5s7b2 


reduced set 


GR hydro - 


h xCFC 


474 ms 


2.8° 


no 




LS180 


M15 


sl5s7b2 


full set 


Newtonian + modified potential 


517 ms 


2.8° 


no 




LS180 


N15 


sl5s7b2 


full set 


Newtonian (purely) 


525 ms 


1.4° 


no 




LS180 


G25 


s25.0 


full set 


GR hydro - 


h xCFC 


440 ms 


1.4° 


no 




LS220 


G27 


s27.0 


full set 


GR hydro - 


h xCFC 


765 ms 


1.4° 


yes 


209 ms 


LS220 



''Defined as the point in time when the average shock radius (rgij) reaches 400 km. 



In other words: At a given time t, x is obtained by di- 
viding the weighted wavelet transform / |x(t, /)| by its 
maximum over all frequencies for that time slice. The 
frequency weighting has been chosen in correspondence 
with the energy spectrum in Equation (I7])I3 We found 
that this weighting and normalization procedure helps to 
reveal the frequency structure of the signal most perspic- 
uously. 

4. EXEMPLARY DISCUSSION OF GW SIGNAL 
FEATURES FOR II.2M0 AND ISM© MODELS 

4.1. Qualitative Description of the GW Signal 

GW amplitudes (both for the matter and neutrino 
signals) for models Gil, G15, M15, N15, and S15 are 
given in the left panels of Figures [JHSl and the evo- 
lution of the signal in frequency space is illustrated 
with the help of wavelet spectra in the right panels of 
these Figures. Qualitatively, our gravitational wave- 
forms share almost all the characteristics of waveforms 
recently obtained from (pseudo-) New tonian simulations 
of core-coll apse supernova explosions ("Marck et alll2009l : 
iMurphv et al. 2009; Yakunin ct al. 2010): We observe a 
quasi-periodic signal during the phase of prompt post- 
shock convection roughly between 10 ms and 50 — 70 ms 
after bounce, which is followed by a more quiescent phase 
of several tens of milliseconds until hot-bubble convec- 
tion and increased SASI activity set in and produce a 
strong stochastic signal component. The stochastic sig- 
nal in Figure [5] is strongest during the first 200 ms af- 
ter the onset of the explosion, and then continues at 
a lower amplitude with proto-neutron star convection 
taking over as the dominant source for high-frequency 
gravitational waves. Model G15 (Figure [3|) also shows 
a "tail" with a steadily increasing wave amplitude in 
the explosion phase, a feature which is due to the ex- 
pansion of a s trong ly prolate shock as recog nized by 
IMurphv et all pOOl and lYakunin et al.l (|20Tol) . In ad- 
dition, anisotropic neutrino emission gives rises to an 
almost monotonically growing amplitude from about 
200 ms onward, which is, in fact, much larger than the 
"tail" in the matter signal. Model Gil, on the other 
hand, exhibits a new feature: Here we observe several 

^ A weighting factor / instead of is used because of the 
usual scale-dependent (i.e. frequency-dependent) normalization of 
the wavelet transform. 



bursts of high-frequency gravitational wave emission in 
the explosion phase with amplitudes comparable to the 
phase of hot-bubble convection - a phenomenon that 
turns out to be connected to the fallback of material 
onto the proto-neut ron star in a rather weak explosion 
(IMiiller et al.ll2ni2bD . 

The wavelet spectra in Figures [2H6] equally reflect the 
evolution through these different phases. The early 
quasi-periodic signal initially produces a peak at ~ 
100 Hz, which is clearly discernible during the first 
^ 100 ms. During the subsequent phase of reduced 
gravitational wave activity, there is no clearly identi- 
fiable narrow emission band, instead broadband low- 
frequency (model Gil until ~ 200 ms) or high-frequency 
(model G15 around ~ 150 ms) noise dominates the spec- 
trum. For models M15, N15, and S15 this intermedi- 
ate phase is not very pronounced in the wavelet spec- 
tra. The subsequent phase of strong SASI activity and 
hot-bubble convection is mostly characterized by a rel- 
atively narrow emission band, which appears between 
100 ms and 200 ms (depending on the simulation) and 
gradually shifts to higher frequencies at later times. Ex- 
cept for model N15, there is also considerable gravi- 
tational wave activity at frequencies above this band. 
This contribution disappears in the later explosion phase 
(> 600 ms), leaving a single, relatively sharply defined 
peak frequency. Note that low-frequency components 
(i.e. asymmetric shock expansion) are poorly reflected in 
the wavelet spectrum due to their small weighting factor 
(see Equation 

In the following, we shall discuss specific features of 
the individual phases of gravitational wave emission in 
our simulations and study the quantitative impact of 
GR on the GW signal in more detail. A cursory glance 
at Figures [2]-|6] already reveals that the spectra contain 
"cleaner" information about the specific signal properties 
of the different models than the stochastically varying 
amplitudes, but we will nonetheless examine the signals 
extensively both in the time and in the frequency do- 
main. 

4.2. Early Quasi-periodic Signal 

With regard to the early quasi-periodic signal, two ma- 
jor questions ought to be addressed, namely whether GR 
affects the typical amplitude and frequency of this signal 
component, and, even more basically, what is the hydro- 
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Figure 2. Left panel: Matter (solid lines) and neutrino (dashed lines) gravitational wave signals for model Gil. Note that the scale for 
the matter signal (left vertical axis) is different from the scale for the neutrino signal (right vertical axis). Right panel: Normalized wavelet 
spectrum x{t> f) of the matter signal for model Gil. The grayscale ranges from white (x = 0) to black (x = 1, maximum value). For the 
definition of xC^i/): see Equation Jsjl. 
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Figure 3. Left panel: Matter (solid linos) and neutrino (dashed lines) gravitational wave signals for model G15. Note that the scale for 
the matter signal (left vertical axis) is different from the scale for the neutrino signal (right vertical axis). Right panel: Normalized wavelet 
spectrum x(t, /) of the matter signal for model G15. 



dynamical origin of the GW signal during this phase? 
Different authors have ascribed the quasi-periodic sig- 
nal during the first several tens of millisec onds directly 
to prompt post-shock con vective overtur n (IMarek et al 



20091: iMurphv et aIll2009D . while others (jYakunin et al 



20101) hav e emphasized t h e con tribution of early SASI 
motions. lYakunin et alj ()2010l ). in particular, argued 
that the quasi-periodic signal is largely due to the decel- 
eration of the infalling matter at an oscillating aspherical 
shock. 



4.2.1. Origin of the Signal 

For our simulations, the question about the origin of 
the "prompt convection" signal can be answered rela- 
tively clearly. The deceleration of matter at the shock 
can be ruled out as major contribution factor with 
the help of an analytic estimate of the expected (mat- 
ter) gravitational wave amplitude shock from a non- 
stationary aspherical shock wave. For' weak SASI oscil- 
lations, we obtain the following formula in terms of the 
pre-shock density pp, the ratio of the post-shock and pre- 
shock densities /3, the power-law index of the pre-shock 
density profile 7, and the multipoles ai of the angle- 



dependent shock position, 
256^3/2^ 



A. 



E2 

20, shock 



5V15c 



Pp - 1) Og [(4 -I- 7)0200 -I- 0200] , 

as derived in Appendix [BJ The multipole coefficients ai 
are defined in terms of the angle-dependent shock posi- 
tion and the ^-th Legendre polynomial as 



2e + i 



rshi0)Pi{O) d cos I 



(10) 



With typical values for the early post-bounce phase 



(Pp 



10^ g cm 



7 



-1.5), the resulting values are 



more than an order of magnitude smaller than the ob- 
served amplitudes on the order of several 10 cm. 

However, convection can be excluded as the direct 
source of the quasi-periodic signal as well, since the 
destabilizing entropy and lepton number gradients are 
erased very quickly, while the initial phase of gravita- 
tional wave emission lasts for several tens of millisec- 
onds. The actual source of the signal can be determined 
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by considering the integrand 
d 

iP = r^ sin 6*— [^^(Scos^ 6-1)- 'dSe sine* cos6'] 4^ 

(11) 

in the quadrupole formula ([T|) for the matter signal in 
order to identify t he main contributions to A^q (cp. 
iMurphv et"aLll2009[ ). By visualizing V (left half of Fig- 
ure [T]), waves or hydrodynamic mass motions responsible 
for gravitational wave emission can be identified fairly 
well. The largest contribution to Afg comes from coher- 
ent stripe-like patterns between 25 km and the shock. 
The underlying flow pattern appears to be dominated 
by propagating wavefronts (and not convective plumes), 
which emerge even more clearly when the partial time 
derivative a = dvr/dt of the velocity field is plotted (right 
panel of Figure [7]). These can be identified as acoustic 
waves by their propagation speed and by the fact that 
the temporal variations 5P and 5p of the pressure and 
the density obey the relation 5P/ P « T5p/ p (where F is 
the adiabatic exponent). The frequency of these waves is 
of the order of 100 Hz, which also accords well with the 
gravitational wave spectrum during this signal phase. In 
our model, prompt convection is thus only indirectly re- 
sponsible for the quasi-periodic signal, either by directly 
initating acoustic waves, or by instigating SASI activity, 
which also involves the propagation of acoustic waves in 
the post-shock region with the SASI frequency. 



4.2.2. Effect of the GR Treatment on the Signal 

The amplitude and frequency of the quasi-periodic sig- 
nal varies appreciably between the models investigated 
here. For the ISM© models with a different treatment 
of gravity (G15, M15, N15), the maximum amplitude 
ranges from 26 cm (G15) to < 2cm in model M15, where 
hardly any gravitational wave activity takes place dur- 
ing this phase. We believe, however, that these differ- 
ences in amplitude may simply stem from a stronger or 
weaker excitation of ^ = 2 shock oscillations and acous- 
tic waves in the different models. The differences are 
therefore not indicative of a systematic effect of the GR 
treatment on the strength of this signal. This conclu- 
sion is also supported by the fact that t he early sig- 
nal is much stronger in model M15LS-2D of IMarek et al.l 
(|2009[) (which only differs from model M15 by a higher 
angular resolution) than in model M15. Overall, the am- 
plitudes in the GR ca se appear to be of similar magni- 
tude as in Newtonian (jKotake et al.ll2007t Murpli v et al. 
2009^ and pseudo-N ewtonian simulations (IMarek et al. 
2009; Ya kunin et al. 2010) with peak amplitudes Afp of 
- 40 cm (Gil) and ~ 30 cm (G15). 

On the other hand, the typical signal frequency is less 
affected by such stochastic variations, and is therefore 
a better indicator for systematic differences, which are 
indeed observed: The gravitational wave spectrum of 
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Figure 6. Left panel: Matter (solid lines) and neutrino (dashed lines) gravitational wave signals for model S15. Note that the scale for 
the matter signal (left vertical axis) is different from the scale for the neutrino signal (right vertical axis). Right panel: Normalized wavelet 
spectrum x(t, /) of the matter signal for model S15. 




Figure 7. The dimensionless integrand ip in the quadrupole for- 
mula |(TJ for the matter signal (left half of figure) and the time 
derivative a = dvr/dt (right half) of the radial velocity field 22 ms 
after bounce for model G15. 
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Figure 8. Gravitational wave energy spectra (matter signal only) 
for models G15 (black), M15 (red) N15 (blue), and S15 (brown) for 
the time interval from 20 ms to 520 ms after bounce. The Fourier 
transform has been carried out without a window function in order 
to retain the high-frequency contribution from the phase of strong 
gravitational wave emission after 400 ms in model G15. In order 
to better differentiate the curves, the spectra have been rescaled 
by a factor of 10-\ lO'^, and 10 for model M15, N15, and S15, 
respectively. 
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Figure 9. Maximum Brunt- Vaisala frequency /p in the con- 
vectively stable region between the proto-neutron star convection 
zone and the gain layer as a function of time for models G15 (solid 
line), M15 (dashed) and N15 (dotted), /p is a rather flat function 
of radius in this region, and the maximum value may therefore be 
taken as an estimate for the "typical" value of the Brunt- Vaisala 
frequency. The plot also shows thepeak frequencies extracted 
from the wavelet spectra (Figures [S] [4] [Sj at intervals of fa 50 ms 
(squares: model G15, circles: M15, triangles: N15). 

the early signal peaks at a somewhat higher frequency 
of sa 100 Hz in the GR case (Gil and G15) compared 
to w 60 Hz in the Newtonian case (N15). For model 
M15, the early signal is rather weak, but there is still a 
peak in the region around « 70 Hz in the wavelet spec- 
trum, which i s in a greement with the result obtained by 
iMarek et all (I2009D fo r their ISM© model with the EoS 
of iLattimer fc Swestvl ()1991[ ). There thus appears to be 
a tendency towards higher frequencies in the GR case. 
This frequency shift is probably the result of a differ- 
ent width and location of the region affected by prompt 
post-shock convection (cp. IMarek eriil. 2009 for this line 
of reasoning in the context of EoS effects on the grav- 
itational wave signal): While prompt convection devel- 
ops in the region between enclosed masses of 0.61Mq 
and O.77Af0 in model G15, the corresponding range in 
model N15 is 0.66Mq . . . 0.83Mq and O.6IM0 . . . 0.71Mq 
in model M 1 5 and in the Lattim er fc Swesty run o f 
IMarek et all (|2009D . Interestingly, IMarek et all (|2009D 
found a similar shift towards higher frequencies with the 
stiffer nuclear equation of state of iHillebrandt fc Wolfl 
(fl985.) . for which the spectrum of the early gravitational 
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wave signal is remarkably similar to that of model G15. 

The early GW signal from model S15 with simplified 
neutrino rates also shows a different frequency structure 
than that of G15. The dominant frequency is initially 
rather high 200 Hz), but after 50 ms post-bounce, 
lower frequencies also appear in the spectrum. Again, the 
width of the layer affected by prompt convection is prob- 
ably a factor responsible for this difference: In S15, there 
are two unstable regions, namely O.77A/0 . . . O.88A/0 and 
0.96Mq . . . 1.O7M0, i.e. the convective region contains a 
significantly larger mass and is located further outside 
than in G15. It is noteworthy that the simplified neu- 
trino rates change the entropy and lepton number pro- 
files in the early post-bounce phase so drastically that 
the dynamics of prompt convection and the GW signal 
are altered significantly. 

4.3. GW Signal from Hot-Bubble Convection and the 

SASI 

4.3.1. Model Comparison - Shift of Characteristic 
Frequencies 

After the early quasi-periodic signal has subsided, a 
phase of relatively weak GW emission ensues, and even 
when hot-bubble convection starts, the GW wave activ- 
ity may not increase immediately. The onset of stronger 
GW emission from hot-bubble convection appears to be 
strongly model-dependent with no clear connection to 
the dynamics. In model Gil, the explosion is already 
underway when the GW amplitude starts to increase sig- 
nificantly 200 ms after bounce, and in the Newtonian 
model N15, there is a similarly long delay. At the other 
end of the scale, we find model S15, where there is no 
gap between the early quasi-periodic signal and the sig- 
nal from hot-bubble convection (which is helped by the 
fact that the gain region already develops at ^ 60 ms in 
this model). 

Due to the stochastic nature of the signal, the ampli- 
tudes of the different models should be compared with 
some caution. All models exhibit amplitudes on the or- 
der of several 10 cm, which is roughly consistent with 
earlier pseudo- Ne wtonian stud i es with sophisticated neu- 
trino transport (iMiiller et all 120041: iMarek et all 120091: 
lYakunin et al.ll2010[ ). Nevertheless, two general trends 
emerge: As exemplified by models Gil and G15, the 
phase of strongest gravitational wave emission begins 
when the shock starts to expand again and lasts. This 
phase lasts for about 200 ms, after which point the 
stochastic signal largely subsides. Furthermore, the grav- 
itational wave amplitude is loosely correlated with the 
mass in the gain region and the convective energy, which 
is why the pessimistic, non-exploding model M15 is char- 
acterized by somewhat smaller amplitudes than G15, a 
mod el with stronger con vection and a larger gain region 
(cp. IMiiller et al.ll2012b[) . For the same reason, model 
S15 with simplified neutrino rates and less favorable 
heating conditions also shows smaller amplitudes than 
model G15 in general. 

As for the early quasi-periodic signal, the gravitational 
wave spectra reveal the effect of general relativity much 
more clearly than the amplitudes. A sizable frequency 
shift in GR compared to the Newtonian and the effective 
potential approximation can be seen both in the wavelet 
spectra in Figures [3]-[5] as well as in the time-integrated 



spectral energy distribution dE/df (see Equation [7]) for 
the first 500 ms after bounce shown in Figure [8] In the 
GR run, the peak is clearly located at a higher frequency 
than in the purely Newtonian case, but the peak fre- 
quency remains somewhat smaller than with the effec- 
tive potential approach. To quantify these differences, 
we compute the median /m of the spectral energy distri- 
bution dE/df, which is implicitly defined by the condi- 
tion 



i.e. fM is the frequency below which half of the total en- 
ergy in gravitational wave is radiated. We obtain values 
of fM = 920 Hz for model G15, 1100 Hz (-1-20%) for 
model M15, and 510 Hz (-44%) for the purely Newto- 
nian model N15. The same ordering is observed in the 
wavelet spectra (Figure [SHS]) throughout the simulation 
once hot-bubble convection starts. As for the early quasi- 
periodic signal, general relativistic effects therefore con- 
siderably affect the gravitational wave spectrum, leading 
to significantly higher frequencies (by ^ 80%) than in the 
purely Newtonian case, while the effective potential ap- 
proach proves to be accurate to within ^ 20%. Again the 
effects of general relativity prove to be of similar mag- 
nitud e as EoS effects (for which we refer to IMarek et al.l 
12001 and emerge as a major factor in determining the 
gravitational wave spectrum. 

Unlike for the early quasi-periodic signal, the effect of 
the simplified neutrino rates is not very pronounced. For 
model S15, we obtain a value of /m = 840 Hz, which is 
somewhat lower than for model G15. At a given time, 
the dominant frequencies are relatively similar for both 
models (FiguresOE]), so the crucial factor for the shift of 
the median frequency must be the onset of the explosion 
in model G15: This leads to enhanced GW emission after 
400 ms and hence a higher weighting factor for late-time 
high-frequency emission than for model S15. 

We note that models Gil and G15 show very similar 
trends in their wavelet spectra, but refer the reader to 
Section [5] for a more thorough discussion of the progeni- 
tor dependence. 

4.3.2. Origin of Frequency Shift in CR 

The dependence of the typical frequency of the gravita- 
tional wave signal from hot-bubble convection and SASI 
activity can be understood by considerin g the domi- 
nant e missi on mechanism during this phase. Marek et al.l 
(|2009f) and iMurphv et"aI1 ()2009f) found that anisotropic 
mass motions in the convectively stable region above the 
proto-neutron star surface actually account for the bulk 
of the signal (although these motions are in turn insti- 
gated by convection and the SASI). While iMurphv et all 
(|2009l ) suggested the deceleration of infalling convective 
plums as a direct so urce for the gravitational wave signal, 
IMarek et al.l (|2009[ ) also mention surface g-mode oscilla- 
tions excited by the downflows as a source. Both mech- 
anisms can work in tandem, and one may speculate that 
the continuous narrow emission bands in Figures [2H5] are 
associated with a time-variable, but relatively stable g- 
mode frequency, whereas the deceleration of plumes with 
random frequency and penetration depth is responsible 
for the more noisy part of the spectrum above this band. 
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In either case, an anisotropic, buoyancy-driven flow is 
responsible for the emission of gravitational waves, and 
the typical angular frequency of the signal is therefore 
approximately given by the buoyancy or Brunt- Vaisala- 
frequency N in the convectively stable regi on between 
the n eutrinosphere and the gain radius fas . Murphy et aH 
120091 explicitly verified for their model). As we shall 
demonstrate, the GR treatment systematically affects 
the buoyancy frequency in a way that accounts for the 
observed frequency differences. 

In the Newtonian case, the buoyan cy frequency is giv en 
by the familiar expression (see, e.g.. lAerts et al.l[2010l ). 



1 a$ 

p dr 



1 dP 

dr 



dp 
dr 



(13) 



where $ is the gravitational potential, p the (rest-mass) 
density, Cg the sound speed, and P the pressure. Equa- 
tion (fT3|) is not applicable in the general relativistic case, 
however, because the underlying equations of hydrody- 
namics arc different. The relativistic expression for the 
Brunt- Vaisala- frequency, derived in Appendix [C] for the 
gauge used in Vertex-CoCoNuT, reads 



dac^ 



dr ph(j)'^ 



1 dP dp{l + e) 



dr 



dr 



(14) 



which contains correction terms involving the lapse func- 
tion a, the conformal factor (j), the specific internal en- 
ergy e and the specific relativistic enthalpy h — l + e/c^ + 

In Figure[9l we compare t he real frequency fp = N/2tt 
(the "plume frequency" of iMurphv et al.l I2009D corre- 
sponding to the buoyancy frequency in the convectively 
stable neutron star surface region for models G15, M15, 
and N15 and also show the peak frequency extracted 
from the wavelet spectrum x at selected points in time. 
We clearly observe the same ordering of the models for 
fp as for the peak frequency, although fp generally over- 
estimates the actual frequency of the spectrum by up to 
30% (as already noticed by Murphv et al. 2009)0 

The different buoyancy frequencies in G15, M15, and 
N15 thus account nicely for the effect of the GR treat- 
ment on the GW spectrum, and the terms responsible for 
the shift of fp can be readily identified: In the Newtonian 
approximation (model N15), the neutron star is less com- 
pact than in GR due to the lack of non-linear strong-field 
effects, and the gravitational acceleration term in d^/dr 
Equation is therefore smaller than the corresponding 
term da(? jdr in GR (or d^/dr in the effective potential 
approximation), fp, and hence the GW peak frequency is 
underestimated. The different compactness of the proto- 
neutron star is most likely also responsible for the dis- 
crepancy between o ur Newtonian results and those of 
IMurphv et al.l (|2009[ ). who obtained even lower GW fre- 
quencies (e.g. only 200 . . . 300 H z at 500 ms): Their use 
of a stiff EoS ()Shen et al.lll998D and of a parametrized 
neutrino cooling and heating scheme, which does not al- 
low for the loss of energy and lepton number from the 

^ The fact that the GW peak frequency is consistently lower than 
fp is consistent with the hypothesis that the dominant frequency 
seen in the GW spectrum is actually that of a surface g-mode 
oscillation, since the buoyancy frequ ency provides an upper bound 
for the g-mode f requency (see, e.g. lKippenhahn fc WeigertI 119901 : 
lAerts et ani20T0l ). 



PNS core, both tend to underestimate the contraction of 
the proto-neutron star. 

On the other hand, the effective potential approach 
(model M15) gives a very good approximation for 
the compactness of the proto-neutron star [d^/dr sa 
dac^ /dr), but still fails to reproduce the correct peak 
frequency because the effects of relativistic kinematics 
are not taken into account. In GR, the correction fac- 
tor ah'^cj)^'^ in Equation (IT4l) reduces fp by about 15% 
to 20%, which largely explains the lower frequencies in 
model G15 compared to M15. In addition, the slightly 
higher neutrino luminosities and mean energies in model 
G15 also correlate with a somewhat more shallow den- 
sity gradient in the neutron star surface region and thus 
reduce fp and the typical gravitational wave frequency 
even a little further. 

These arguments indicate that gravitational wave spec- 
tra from the later (> 150 ms) post-bounce phase with an 
accuracy better than ^ 20% in frequency space can only 
be obtained within the framework of general relativis- 
tic hydrodynamics, since both the Newtonian and the 
pseudo-Newtonian approximation suffer from an intrin- 
sic accuracy limit. 

4.3.3. Relation between PNS Properties and the 
Characteristic GW Frequency 

Could the characteristic GW frequency, which emerges 
clearly from the wavelet spectra in Figures [3H5] provide 
clues about properties of the proto-neutron star, such as 
its mass, compactness, or surface temperature? A simple 
analytic estimate for fp based on Equation ([T4|) can shed 
light on this question. 

The metric functions a and in Equation (jl4p can be 
approximated in terms of the proto-neutron star mass M 
and radius i? as a « Ina « « 1 — GM/{Rc^) at the 
proto-neutron star surface, and the pressure and density 
gradients can be found by assuming a roughly isother- 
mal stratification (with temperature T) in the convec- 
tively stable neutron star surface layer. Furthermore, an 
ideal gas equation of state can be used as non-relativisticQ 
baryons dominate the pressure in the relevant region. 
The gradient terms in Equation ()14p can then be imme- 
diately obtained in terms of T and the neutron mass m„ 
from the equation of hydrostatic equilibrium. 



dP _ d f pkT 
dr dr 



9 In a 

dr 



(15) 



Employing the relation — TkT/mn for the speed of 
sound (r being the adiabatic index), we then arrive arrive 
at the following expression for fp, 



1 GM /(r-l)m^ 



N 



27r 27r 



ThT 



1 



(16) 



where F is the adiabatic index in the proto-neutron star 
surface region, and to„ is the neutron mass. The mean 
energy {Ep^ ) of electron antineutrinos may be used as a 
proxy for the temperature (with additional redshift cor- 
rections, i.e. (-Ep^) ~ 3.15lQ;(i?)T), and one thus obtains 
a fairly accurate formula for the evolution of the domi- 

^ This implies e <C 1. 
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Figure 10. Comparison of the wavelet spectra (identical to those in Figures [2] and [S] respectively) of model Gil (lcft)and G15 (right) and 
the prediction of Equation I I17I1 for the typical gravitational wave frequency in terms of the proto-neutron star mass and radius and the 
mean energy of electron antineutrinos measured by an observer at infinity. Equation ljl7|l describes the evolution of the typical frequency 
fairly well and only overestimates /peak somewhat at late times. 



nant GW frequency /peak: 



peak 



1 GM 



1.1 



1 - 



gmv 



(17) 



For our two exemplary models Gil and G15, Figure [TUl 
shows that Equation ([T71) is in fairly good agreement 
with the measured frequencies. The critical parameters 
regulating /peak are thus i) the surface gravity, ii) the 
surface temperature, and iii) the compactness parame- 
ter GM/Rc^ of the proto-neutron star. Since /peak also 
depends on the thermal properties of the PNS surface 
layer, gravitational waves from the accretion phase prob- 
ably cannot provide an unambiguous probe for the bulk 
properties of the proto-neutron star (mass, radius). 

Equation (flT)) illustrates that /peak is sensitive to fac- 
tors that affect the contraction and thermal evolution of 
the proto-neutron star, such as the EoS and the neu- 
trino treatment. The different neutrino rates in model 
S15 are not critical in this respect since they affect the 
proto-neutron star surface temperature only on a very 
moderate level (see paper II). More radical approxima- 
tions in the neutrino treatment could potentially have 
a sizable impact, however, provided that they change 
the proto-neutron star surface temperature and the neu- 
trino mean energies considerably. The neutrino treat- 
ment (multi-group variable Eddington factor transport 
vs. parametrized heating and cooling) along with the 
different equation of state may also partially account for 
the d ifferences of ou r Newt onian model N15 to the mod- 
els of iMurphv et al.l (|2009l ). For the dependence of /peak 
on the progenitor, we refer the reader to Section [5] 

4.4. Explosion Phase ~ Proto-Neutron Star Convection 
and Late-Time Bursts 



As established by IMurphv et all ()2009[ ) and 



lYakunin et all (|2010D . the emission of high-frequency 
gravitational waves subsides to a reduced level once the 
shock accelerates outwards because the excitation of 
oscillations in the proto-neutron star surface by violent 
convective motions in the hot-bubble region largely 
ceases. This is well reflected in the GW signals of 
models Gil and G15 (Figures [2l [S]) in the initial phase 
of the explosion, as the GW amplitude drops noticeably 
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Figure 11. Energy Bgw radiated in gravitational waves as func- 
tion of time for the explosion models Gil (black solid line) and 
G15 (red dashed line) for the entire duration of the simulation. 
For model G15, most of the energy is radiated around the onset of 
the explosion between 400 ms and 600 ms after bounce. By con- 
trast, late-time GW bursts carry a sizable fraction of the radiated 
energy in the case of model Gil. 



after 400 ms and 600 ms, respectively. However, the 
surface g-mode oscillations remain the dominant source 
of high-frequency gravitational waves during this phase: 
The wavelet spectra clearly show that the late-time 
signal comes from the same emission band as the signal 
from the accretion phase (which shifts continuously to 
higher frequencies as the neutron star contracts). An 
analysis of the integrand ip in the quadrupole formula ([1]) 
as in Section 14.2.11 also confirms that aspherical motions 
in the proto-neutron star surface region are still the 
dominant source of gravitational waves. Proto-neutron 
star convection now provides an excitation mechanism 
for the oscillations ([Miiller k Jankal llQQTt iMarek et"all 



[2U09I: livfarohv et al.l [20091 iMiiller et al.il2012c[ ). but due 

the subsonic character (with Mach numbers < 0.05) of 
the convective motions, the GW amplitude remains 
rather small. 

However, model Gil contradicts the established pic- 
ture of subsiding GW emission during the explosion. For 
this progenitor, we observe several "bursts" of stronger 
gravitational wave emission later in the explosion phase 
on at least three occasions (630 ms, 790 ms, and 820 ms). 
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Figure 12. Snapshots of the entropy s in the region around the proto-neutron star at post-bounce times of 812.3 ms, 816.5 ms, 819.7 ms, 
and 821 ms, depicting the formation of a new downfiow by material falling through a hot bubble of neutrino-heated material. To guide the 
eye, the dent in the high-entropy bubble, which eventually develops into the new downfiow is highlighted by a white arrow in each panel. 
The values of the entropy range from Ofeft/nucleon (black) to 35A;(,/nucleon (yellow). Due to its high infall velocity, the downfiow overshoots 
far into the convectively stable proto-neutron star surface layer, exciting violent oscillations. The GW burst occurring at this time can be 
seen in Figures [2l and [T3l 
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Figure 13. Correlation of late-time bursts in the GW signal and 
the neutrino signal; The top panel shows the GW amplitude 
as a function of time for model Gil between 600 ms and 920 ms 
after bounce. The middle panel displays the neutrino luminosity 
of Ue (black solid line), Pe (red, dashed), and i^^/t- (blue, dash- 
dotted) - defined as the neutrino energy flux integrated over all 
emission directions as measured at an observer radius of 400 km. 
The bottom panel shows the direction-averaged mean energy of 
neutrinos at the same observer radius. 

during which the amplitude can become comparable to 
that coming from the phase of hot-bubble convection. 
About 40% of the GW energy is emitted later than 
600 ms after bounce (Figure [TT|) . which is in stark con- 
trast to the more vigorous explosion in model G15 with 
faster shock expansion. These bursts occur when matter 
falls back onto the proto-neutron star through a newly 
developing downfiow, which excites g-mode oscillations 
in the proto-neutron star surface region (Figure fT2|) . The 
formation of new downflows in the explosion phase is 




Figure 14. Gravitational wave energy spectra (matter signal 
only) for the relativistic models G8.1 (orange), G9.6 (magenta), 
Gil (red), G15 (black), G25 (green), and G27 (blue), for the entire 
duration of the simulation. Note that some of the spectra have been 
rescaled by a factor given along with the model designation. The 
Fourier transform has been carried out without a window function 
in order to retain the high-frequency contributions towards the end 
of some simulations. The spectra reflect larger amplitudes of the 
early quasi-periodic signal (around ~ 100 Hz) in models Gil, G25, 
and G27. For model G9.6, there is almost no high-frequency sig- 
nal from hot-bubble convection. For Gil, the sharp peak slightly 
above 1 kHz emerges in the spectrum due to late-time GW bursts. 



a consequence of the small explosion energy of model 
Gil (see paper II), and one could speculate that low- 
energy fallback supernovae may generally reveal them- 
selves through such multiple GW burst episodes. 

Several interesting features of these late-time bursts 
should be noted: Between the few instances where new 
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accretion downflow develop, aspherical mass motions 
above the proto-neutron star surface are apparently not 
strong enough to produce significant noise in the GW sig- 
nal, and the g-mode frequency therefore emerges much 
more cleanly from the spectrum during the late time GW 
bursts than in the accretion phase. The narrow-band 
character of the spectrum (see Figure[2]and the red curve 
in Figure fT4)) is potentially helpful both for the detection 
of the GW signal as well as for the interpretation of the 
data (e.g. by allowing for a clearer discrimination of dif- 
ferent equations of state on the basis of the g-niodc fre- 
quency). Moreover, the bursts directly coincide with an 
enhancement of the electron neutrino and antineutrino 
luminosities and mean energies as shown in Figure 1131 
When the newly-formed downflows bring fresh material 
into the cooling region, the luminosity of i^e and Pe in- 
creases abruptly by up to ^ 20%, and the mean energy 
of the emitted neutrinos jumps by up to ~ 1 MeV. Such 
correlations between the neutrino and GW signals prob- 
ably merit further investigation, but a deeper analysis 
taking into account both the observer position and the 
detailed anis otropics in the neut rino radiation field in 
the manner of^Mii ller et al.l (j2012d ) will not be attempted 
here. The anisotropic neutrino emission will be discussed 
more thoroughly in a subsequent publication. 

4.5. The Tails in the Matter and Neutrino Signals 

At late times, the gravitational wave signals of our 
models Gil and G15 exhibit the same low-frequency fea- 
tures that have been observed in previous studies: The 
matter signal develops an offset or "tail" which has been 
ascribed to aspherica l shoc k propagation (jMurphv et al.l 
I2009t lYakunin et"all l2010f ). with prolate/oblate explo- 
sions leading to positive/negative amplitudes. This offset 
is of the order of 10— 15 cm for model G15 (which is con- 
sistent with its prolate explosion geometry) , while only a 
small offset of ~ 3 cm develops for model Gil at rather 
late times due to the small shock deformation. These 
amplitudes are consistent with a direct numerical evalu- 
ation of Equation (jB3[F I. which confirms that aspherical 
shock propagation is indeed responsible for the offset in 
the signal. 

However, as in lYakunin et al.l (|2010f ) the "tail" sig- 
nal is much weaker than the similar tail-like signal from 
the anisotropic emission of neutrinos. Models Gil and 
G15 both develop long-lived polar downflows that sub- 
sist well into the explosion phase and lead to a sustained 
quadrupolar emission anisotropy over several hundreds 
of milliseconds; in model G15 the running average of the 
anisotropy parameter a^, (Equation [5]) over 50 ms can 
become as large as 0.02. This is in striking contrast to 
model M15 with many alternating episodes of enhanced 
emission from the polar and the equatorial region. It 
is very likely, however, that stochastic model variations 
affect the neutrino-generated gravitational wave signal 
considerably; model M15, e. g. , is d istinguished from the 
L&S model of iMarek et all (I2009D o nly by a different 
angular resolution, vet iMarek et ahl (|2009t l obtained a 

* Note that Equations l|9IB6| l cannot be used to evaluate the 
gravitational wave amplitude due to aspherical shock propagation 
for model G15 because the Legendre coefficient ai of the angle- 
dependent shock position r{9) is of th e sam e order as ao, and the 
assumptions used to derive Equations l)9IB6|l therefore break down. 



large negative signal amplitude due to the predominant 
emission of neutrinos from the equatorial region. How- 
ever, even in model M15, the typical amplitude of the 
neutrino-generated signal is larger than the typical sig- 
nal from aspherical shock propagation, which will prob- 
ably always remain a subdominant contribution to the 
low-frequency part of the spectrum. 

5. PROGENITOR DEPENDENCE OF THE GW 
SIGNAL 

After having discussed the different phases of GW 
emission for models Gil and G15, we now turn to the 
progenitor dependence of the GW signal. An overview 
over the waveforms for models G8.1-G27 is given in Fig- 
ure 1151 and time-integrated energy spectra are shown in 
Figure 1141 Table [5] summarizes important characteristic 
numbers for the waveforms and spectra. 

5.1. Progenitor Dependence of Waveforms 

Figure [15] and Table [2] illustrate a few general trends 
in the GW signals: GW activity tends to increase for 
models with higher accretion rate M during the phase 
of convective and SASI activity. Table [2] shows that ex- 
cept for model G25 there is a strong correlation between 
the accretion rate Mgo at a post-bounce time of 90 ms 
(typically close to the onset of strong convection/SASI) 
and the total gravitational wave energy -Egw well as 
the maximum amplitude ^fo max- This correlation is rea- 
sonable, since the typical GW amplitude depends on the 
mass in the gain and cooling regions participating in vio- 
lent aspherical motions and g-mode activity, respectively. 
Both of these masses in turn regulated by the accretion 
rate. Naturally, Mgo should not be understood as a sin- 
gle unambiguous parameter regulating the GW emission; 
it is rather the overall evolution of M that is relevant. 

Because of vastly different accretion rates, models G9.6 
and G27 constitute two extreme cases on the scale of 
GW emission: In model G9.6, the accretion rate drops 
so quickly that there is almost no signal from hot-bubble 
convection, and only the early quasi-periodic signal re- 
mains. Moreover, although there is some convective over- 
turn in the ejecta in this model, the rather rapid ex- 
plosion largely precludes the excitation of PNS surface 
g- modes by downflows or accretion downdrafts. The 
matter affected by overturn is blown out rather quickly 
ahead of the developing neutrino-driven wind (cp. top 
right panel of Figure [T]) . 

By contrast, a lot of mass is involved in strong SASI os- 
cillations prior to the onset of the explosion in model G27 
(bottom right panel of Figure [ij. Moreover, the shock 
does not expand as rapidly afterwards, thus allowing vi- 
olent aspherical motions to continue for several hundreds 
of milliseconds beyond the onset of the explosion. Con- 
sequently, the energy emitted into gravitational waves is 
as high as 2 X 10^^ erg for this moodel. 

While the models G8.1, Gil, and G15 fit nicely into 
a continuum between these extreme cases, model G25 
shows less GW activity than model G27 although it is 
characterized by even higher mass accretion rates as well 
as by continuing SASI activity (bottom left panel of Fig- 
ure [T]) . This reversal of the aforementioned trend is re- 
lated to fact that G25 does not develop an explosion 
within the simulation time. GW emission is typically 
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Table 2 

Progenitor dependence of the GW signal 



Model 


M90" 


(10-*^ erg) 


AH2 c 
^20, max 

(cm) 


(Hz) 


(Hz) 


/40()' 
(Hz) 


(Hz) 


MpNs'^ 

(Mq) 


^cxpl 

(ms) 


G8.1 


0.38 


0.36 


45 


540 


600 




730 


1.37 


177 


G9.6 


0.09 


0.0032 


13 








130 


1.36 


125 


Gil 


0.57 


2.7 


52 


570 


580 


680 


1100 


1.35 


213 


G15 


1.08 


4.8 


70 


520 


580 


710 


950 


1.58 


569 


G25 


1.80 


1.2 


77 


660 


750 


900 


620 


2.01 




G27 


1.44 


20 


253 


580 


700 


860 


680 


1.77 


209 



"Mass accretion rate 90 ms after bounce. 
''Energy of radiated gravitational waves. 
"^Maximum GW amplitude of the matter signal. 
'^Dominant emission frequency 250 ms after bounce. 
Dominant emission frequency 300 ms after bounce. 
'Dominant emission frequency 400 ms after bounce. 
^Median frequency of the GW energy spectrum. 
''Baryonic mass of proto-neutron star by the end of the simulation. 

'Explosion time, defined by the instance when the average shock radius (rah) reaches 400 km. 



strongest around the onset of the explosion (G8.1, Gil, 
G15, and G27), when the mass in the gain region in- 
creases considerably and aspherical motions in the post- 
shock region become most violent (i.e. reach the highest 
velocities) before being quenched again during the later 
phases of shock expansion. This phase is lacking in model 
G25, and the strong retraction of the shock instead re- 
duces the mass in the gain region so that the GW am- 
plitude is strongly attenuated during the later evolution. 
Moreover, the relatively pure £ = 1 SASI mode seen in 
model G25 may not be very efficient in exciting £ = 2 
perturbations in the flow to which the GW amplitude is 
sensitive. 

The overall trends in the progenitor dependence sug- 
gested by the six models G8.1-G27 can thus be summa- 
rized as follows: GW emission tends to become stronger 
with increasing mass accretion rate. In addition, both a 
strong shock retraction and a very rapid explosion can 
quench GW activity. 



5.2. Progenitor Dependence of GW Spectra 

To characterize the spectral properties of models G8.1- 
G27, we show time-integrated energy spectra in Fig- 
ure 1141 foi' which we give the median frequency /m in 
Table [5] Moreover, we list the dominant emission fre- 
quencies /250, /30O: and /400 around 250 ms, 300 ms, and 
400 ms after bounce to show the time-frequency evolu- 
tion of the GW signals. 

For the early quasi-periodic signal, which emerges as 
a low-frequency peak in the curves of Figure 1141 we find 
no significant progenitor dependence. This signal com- 
ponent always peaks around 100 Hz in our models. The 
later signal from convection and the SASI shows some- 
what larger variations. The typical emission frequencies 
/250, /300, and /400 vary by 30%, and tend to be high- 
est for model G25. Considering that we cover a very 
wide range in PNS masses from ~ 1.36M0 (G9.6, Gil) 
to ^ 2.OAf0 (G25), this variation may seem astonishingly 
small. This can be understood as a partial cancellation 
of terms in Equation (|17p for the dominant emission fre- 



quency: 

fpcak 



1 GM 
2^~R2- 



1.1 



(18) 



While the surface gravity GM/R^ is systematically larger 
for more massive and compact neutron stars, the proto- 
neutron star surface temperature is also higher, and the 
general relativistic correction term (1 — GM/{Rc'^))'^ is 
smaller. The remaining net change is typically moderate, 
as illustrated by Gil and G25 as rather extreme exam- 
ples in Figure [16) Although the surface gravity is higher 
by up to ^ 90% in model G25, /peak is only lower only 
by ^ 35% in model Gil. Among the less massive proto- 
neutron stars (models G8.1, Gil, G15), the differences 
are far less pronounced. The vari ation between progen - 
itors is similar to that found by iMurphv et al.l ()2009D , 
although the absolute values for the frequencies are very 
different (since Murphy et al. 2009 lacked self-consistent 
neutrino transport and GR in their simulations, and used 
different progenitors as well as a different EoS). 

When we consider time-integrated GW energy spectra 
(Figure [131), the relation between PNS and GW prop- 
erties becomes more complicated, however. We observe 
the highest median frequency of /m — 1100 Hz for model 
Gil, while we find relatively low values between 600 Hz 
and 700 Hz for G25 and G27, although these two models 
have higher values for /250, /300, and /4oo- This rever- 
sal occurs because much of the gravitational wave emis- 
sion in Gil happens after the onset of the explosion in 
late-time bursts (see Section |4^ . GW waves with rela- 
tively high frequencies therefore contribute more strongly 
to the spectrum than in models such as G25 and G27. 
These late-time bursts are distinctly visible as a peak 
around 1100 Hz for Gil in Figure 111 For model G15, 
the situation is somewhat similar, because the bulk of 
GW emission occurs around 500 ms due to the rather 
late explosion. Model G9.6 constitutes another extreme 
because the high-frequency component of the spectrum 
is largely absent, resulting in a low median frequency of 
fM = 130 Hz. 

Our results demonstrate that time-integrated spectra 
are the result of a somewhat complicated interplay be- 
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Figure 15. Matter and neutrino gravitational wave amplitudes Ajg from GR simulations of the six progenitors considered in this paper. 
The matter and neutrino signals are shown as solid black and dashed red curves, respectively. For models G15 and G27, the gravitational 
wave signal from the anisotropic emission of neutrinos has been rescaled by a factor of 0.2. For the exploding models, the onset of the 
explosion (defined as the time when the average shock radius reaches 400 km) is indicated by a dotted blue line. 
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Figure 16. Comparison of the terms contrib utin g to tfie gravita- 
tional wave peak frequency /peak in Equation I I17I1 for models Gil 
(solid lines) and G25 (dashed lines). The proto- neutron star surface 
temperature T and the temperature (i?^^/3.151) corresponding to 
the mean electron antineutrino energy (assuming a Fermi distri- 
bution with vanishing chemical potential) are shown in the middle 
panel. The estimated neutron star surface gravity GM/B? and the 
the relativistic correction factor (1 — GM/Rx?) in Equation l|17ll 
are shown in the bottom panel, and the resulting estimate for the 
peak frequency is depicted in the top panel. Note that we define 
the neutron star surface by a fiducial density of 10^^^ g cm""^. 

tween the time of strongest GW emission and the time- 
dependence of the dominant emission frequency. This 
impHes that time-integrated spectra are even more dif- 
ficul t to r elate to PNS properties than /peak (cp. Sec- 
tion [4331). 

6. SUMMARY AND CONCLUSIONS 

Based on sever al recent two -dimensional relativistic 
explosion models (jMiiller et al.l [2012b) for non-rotating 
progenitors between S.IMq and 21Mq, we studied the 
GW signal of core-collapse supernovae from the early 
post-bounce phase well into the explosion phase. We 
not only provided gravitational waveforms from GR hy- 
drodynamics simulations with sophisticated multi-group 
neutrino transport for the first time, but also analyzed 
the impact of the GR treatment (GR vs. effective poten- 
tial approach vs. Newtonian approximation) and of the 
neutrino physics input on the GW signal by means of 
three complementary (non-exploding) simulations of the 
15Mq progenitor. 

In all phases, we found that GR has a sizable im- 
pact on the GW spectrum. For purely Newtonian mod- 
els, we found that the typical GW frequency during 
the phase of hot-bubble convection is severely under- 
estimated (by 40% . . . 50%) compared to the GR case, 
while it tends to be overestimated by ~ 20% in the ef- 
fective potential approach. We determined that this is 
the results of systematic differences of the buoyancy fre- 
quency in the proto-neutron star surface region (which 
approximates the typical GW frequency as pointed out 
by .Murphv et al.. ,2009.) . In the Newtonian approxima- 



tion, the smaller compactness and surface gravity of the 
neutron star lead to a lower buoyancy frequency com- 
pared to GR, while the effective relativistic potential ap- 
proach fails to correctly reproduce the frequency of oscil- 
lations in the proto-neutron star surface because of miss- 
ing relativistic correction terms in the equations of hy- 
drodynamics. Similarly large differences as for the phase 
of hot-bubble convection were also seen for the quasi- 
periodic signal produced by prompt convection and early 
SASI activity. 

Apart from systematic GR effects, we observed novel 
features in the GW signal during the explosion phase for 
the 11.2Mo progenitor, a model that is likely to repre- 
sents a fallback supernova. Here the fallback of mate- 
rial onto the proto-neutron star through newly-forming 
downflows in the explosion phase leads to a strong ex- 
citation of surface g-mode oscillations, which produce 
burst-like GW signals that are sharply defined in fre- 
quency space. These late-time bursts are correlated with 
small bursts of enhanced Ve and emiss ion. Contrary 
to wh at more energetic exp losion models ([Murphv et al.l 
120091 : lYakunin et all 120101 ) might suggest, strong GW 
emission in the high-frequency band late in the explosion 
phase appears to be possible for underenergetic super- 
novae. In all other respects, however, our models qual- 
itatively confirm the four-phase picture of gravitational 
emission with an early quasi-periodic signal, a quiescent 
phase of several tens of milliseconds, a strong stochastic 
GW signal lasting until some fraction of a second after 
the onset of the explosion, and a low- frequency tail-signal 
from anisotropic neutrino emission and shock expansion. 

Our findings clearly demonstrate that a proper treat- 
ment of GR is indispensable for obtaining accurate pre- 
dictions of gravitational waveforms and spectra. In par- 
ticular, GR simulations will be required to link the GW 
signal to the properties of the nuclear equation of state, 
or of the progenitor. However, as a cautionary remark, 
we emphasize that the inclusion of GR effects at the ex- 
pense of detailed multi-group neutrino transport may po- 
tentially yield little gain in our opinion. Although the 
simplifcations of the neutrino rates considered as test 
case in this paper cause a significant difference only for 
the early quasi-periodic signal, more radical approxima- 
tions in the neutrino transport sector may significantly 
affect the contraction, the compactness, and the surface 
stratification of the neutron star, and might thus severely 
alter the GW spectra. 

Our simulations also revealed trends for the progeni- 
tor dependence of the GW signal. We found a variation 
of the maximum signal amplitude by more than an or- 
der of magnitude between our two extreme cases with 
9.6AfQ and 27Mq. In general, high mass accretion rates 
turn out to be conducive to strong GW activity provided 
that an explosion still develops. The time-dependent 
dominant emission frequency varies by ~ 30% with a 
tendency towards higher values for more massive proto- 
neutron stars. The explosion dynamics plays a major role 
for determining the time-integrated spectrum, which is 
therefore more difficult to relate to PNS properties. 
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APPENDIX 

A. GRAVITATIONAL WAVE EXTRACTION: THE MATTER SIGNAL 

Although the standard Newtonian quadrupole formula for the matter signal and its time-integrated varieties (Finn) 
[l989; .Finn fc Evans. .1990; .Blanchet et al...l99M as used in the CoCoNuT code so far (Dimmelmeier et al. .2002bl 
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l2008f) has proved reasonably accurate (i.e. on a level of 10% . . . 20%) even for relatively strong g r avitational fields i n 
comparison with more sophisticated direct wave extraction methods (jShibata fc Sekiguchil 120031 : iNagar et al.l I2007D , 



- ( «i * j,4~ij ) ■ (Al) 



it suffers from several drawbacks: Its application in the context of general relativistic hydrodynamics simulation is 
beset with a number of ambiguities concerning the proper generalization of the source density p. Moreover, the 
transformation of the second time derivative of the mass quadrupole moment to a numerically more tractable form 
(time-integrated quadrupole formula, stress formula) that does not introduce excessive numerical noise becomes less 
straightforward than in the Newtonian limit. These problems can gi ve rise to certain patholog ies in the computed 
gravitational wave signals such as spurious offsets (cp. Appendix A in iDimmelmeier et al.|[2002b[ ). most notably (but 
not exclusively) with the Newtonian stress formula. 

Evolving the full space-time metric and using direct extraction would be an obvious solution to the problem, but 
this would not only require that we sacrifice the stable xCFC formulation of the metric, but introduce other difficulties 
as well: Direct extraction methods can be expected to be very vulnerable to numerical noise as the gravitational wave 
signal from core-collapse supernovae is weak in the sense that it carries away only a minuscule frac tion of the total 
mass of the system. Post-Newtonian generalizations of the quadrupole formula (|Blanchet et al.|[l990[ ) provide another 
possible workaround, but tend to be considerably more cumbersome than the standard quadrupole formula. For the 
extraction of gravitational waves from our core-collapse supernova simulations, we prefer a somewhat simpler solution 
that still meets the following requirements: 

• No higher time derivatives are required for extracting the wave signal. 

• The gravitational wave amplitude vanishes for a stationary flow. 

• The generalization of the quadrupole formula should (and need only) be valid for the conditions typically en- 
countered in core-collapse supernova simulations. 

Due to the restriction to the supernova context, several approximations can be made: Concerning the matter sources, 
we assume all velocities to be sufficiently small to disregard retardation effects. Concerning the space-time metric, we 
consider only small perturbations on the background of a CFC metric 17^1,, which is given in terms of the lapse function 
a, the shift vector the conformal factor and the flat-space three metric as 

- A/3* 

We also assume that the space-time is almost stationary with vanishing shift, 

t»0, f ,-.„^ ,A2, 

Furthermore, we observe empirically that for the typical field strength reached in the proto- neutron star environment, 
the lapse function and the conformal factor obey the relation 

a<?!)2 w 1, (A3) 

which is a direct consequence of the fact that the field strength is moderate and that the Eulerian energy density E 
is the dominating source term in both the non-linear Poisson equations for a4> and (j) (cp. paper I). 

With these approximations for the space-time metric, it can easily be verified that the harmonic gauge con dition 
d/dx'^{-J^g^'''^) = is satisfied. We can therefore use the relaxed field equations (|Straumaiinll2004l : iPati fc Willll2000) 
to compute gravitational amplitudes perturbatively on the background of the CFC metric. The relaxed field equations 
in geometrized units read 

= 167rG(-5) (T'^'^ + r^l) + V{rn, (A4) 

with g — det((7^y), g^^^ = {—g)g^^- Here, T'^'^ and t^l are the stress-energy tensor of the matter and the Landau- 
Lifschitz pseudo-tensor for the gravitational field, respectively, and I? is a non-linear differential operator. 

The relaxed field equations are amenable to an iterative solution by means of the Green's function method, i.e. 
solutions for /i"^ can be obtained (jStraumannll2004 ) by successively computing 

..■-ftx) = 4 ''-'--''■"'' dV. (A5, 
J |x-x'| 

and updating the source r^" = {—g){T^'' + t^l) + 2'(3^'')/(167rG) in the process. In general, this procedure is beset 
with a number of difficulties (such as the treatment of non-compact sources and the occurrence of divergences in 
higher-order approximations), but it can readily be used to obtain an estimate for the gravitational wave amplitude 
for the matter signal in our case: As in the usual derivation of the Newtonian quadrupole formula, we assume tha t the 
source is small compared to the wavelength and that retardation effects can be neglected, in which case Equation (|A5I) 
can be written a r 

M'''(t,x)«— / T^''(<' =t- IxLxOd^x'. (A6) 
|x| J 

^ We emphasize that unhke Equation IIASI l. Equation IIA6I 1 can , ■, ■ ■■ j-ji- 

, 1 J i i ii .J ,^ ! , to anisotropic neutrino emission, which must include retardation 

no lonser be used to compute the gravitational wave signal due at- t ■ t- -li. ti. 1 i iu 

° t- t> b enects as neutrinos propagate with the same velocity as the gravi- 

tational waves themselves. 
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By projecting h'^'^ onto the subspace of symmetric transverse-traceless (STT) tensors, one obtai ns the angle-dependent 
gravitational wave amphtude, which can then be decomposed into pure-spin tensor harmonics (|Mathewslll962l : fZerillil 
Il970t iThorndlTgSOl) of de gree £ = 2. These are already completely determined by the purely spatial components h^^ . 

At this point, the complicated metric-dependent terms in r'^'^ are the major obstacle for formulating a simple integral 
expression for the gravitational wave amplitude. However, using the fact that r^'' is divergence- free {di^r'^'^ — 0), the 
integral in Equation ()A6P can be transformed into an integral over t'^^ or t^'^ at the cost of introducing time derivatives: 



/i*J'(i,x) 



X dt ' ^ 



J xV"J(t- |x|,x')d3.x' = |||^ J x'x^T"°{t-\x\,K')d^x'. 



(A7) 



The formula for ft,'-' containing the time derivative of t^^ is particularly useful because it does not contain purely 
gravitational contributions due to our assumptions for the metric functions, and can therefore be expressed in terms 
of the covariant Eulerian three-momentum density 5^ = phW'^Vj as 



r"3 



(A8) 



which is exactly the conserved quantity in our formulation of the equations of hydrodynamics (jBanvuls et al.iri997D . 
Plugging in the expression for d^/jSj/dt (cp. also Equation ICip . we therefore obtain 



I dt 



(t-|x|,x')d^x' = 



hyd 



dx^ 



dx^ 



(i-|x|,x')d3a;'. 

(A9) 

Note that the time derivative of -J^Sj includes only the (hydrodynamical) changes due to advection, pressure gradients 
and gravitational source terms, but no neutrino momentum source term ^J^Sjy. 



\ dt 



hyd 



dt 



(AlO) 



As ^/^Sj^u is exactly balanced by a source term for the neutrino energy- momentum tensor, momentum exchange 
between the matt er a nd the neutrinos does not contribute as a source for gravitational waves. 

While formula (jA9|) for h"^^ is already amenable to a numerical treatment (no time derivatives are required), it can 
be cast into a form corresponding more closely to the Newtonian stress formula. Using partial integration, the term 
containing the flux divergence can be eliminated again, 



/i*J(x) 



4 



8x3 



+ {Sjv' + PS)) d- 



^x' 



(All) 



where we have suppressed the arguments t and t' for the sake of simplicity. Upon extracting the symmetric transverse- 
tracefree (STT) part of the right-hand side, the pressure term drops out. 



4 



dxi 



STT 



The resemblance to the Newtonian stress formula, (jNakamura fc Ooharalll989l : iBlanchet et al.|[l99ClD 



(A12) 



(A13) 



STT 



is evident. The decomposition of the gravitational radiation field given by Equation (jA12[) into pure-spin harmonics 
proceeds exactly as in the Newtonian case. In the case of axisymmetry, only the "electric" quadrupole with to = is 
present, and its amplitude Afg can be expressed in spherical polar coordinates as. 



A. 



E2 _ 



32^3/2 



J a<j)^r^ sine [S',.^ (Scos^e* 



1) +Sev" (2-3cos2 



6rSrV sin 6* cos 9 



-Sr,giavr (3cos^ 9 - l) 



3S0 



smWcosf 



d9dr, 



where 5* designates the gravitational source term in the momentum equation: 

2 dx' ■ 



-■ijgrav 



(A14) 



Alternatively, one can retain the time-derivative d^/^Sj /dt in Equation (jA9l) and formulate a modified version of the 
time-integrated quadrupole formula, but the neutrino momentum source term needs to be handled with care in this 
case. In order to avoid double counting, these source terms must be subtracted from the time derivative of y/jSj, and 
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we therefore obtain 



A 



E2 



327r3/2 Q 



yiS dt 

32^3/2 



15 



j (t)^r^ sin e [Sr (3 cos^ 9-1)+ ir^'^Se sin 9 cos ( 
Sr,u^ (3 cos^ 6* — l) + 3r~2S'e,i/ sin 9 cos ^ 



■^r"^ sin 



(A15) 



Since our niodified stress formula has been obtained through exact transformations of Equation (jAlip , which contains 
the time derivative of Sj under the integral, it does not produce any unphysical offset in the gravitational wave signal for 
stationary configurations with dSj/dt = in the limit of infinite spatial resolution^ It also avoids the contamination 
of the matter signal by neutrino momentum source terms since it allows for a nice separation of the source density for 
the wave equation into terms due to the hydrodynamic momentum flux and gravitational source terms. However, the 
contribution from neutrino stresses is small in practice, and may even be neglected in Equation (|Alip with little loss 
of accuracy. 

B. GRAVITATIONAL WAVE SIGNAL FROM ASPHERICAL SHOCK PROPAGATION 

As mentioned in Section 21 it is possible to formulate a simple analytic expression for the contribution of asymmetric 
shock propagation to the gravitational wave signal in the limit of a relatively small deformation of the shock. To this 
end, we perform the integral in the quadrupole formula (jA15l) over a narrow region around the shock between ra{9) 
and rb{9) in radius such as to neglect slower, secular variations in the pre- and post-shock conditions. For the sake 
of simplicity, we disregard relativistic kinematics and strong-field effects, setting a — (f) — 1 and Si = pvi for the 
momentum density. As we confine ourselves to a weakly deformed shock, we also neglect the non-radial component 
V0 of the post-shock velocitjQ and obtain the contribution of the shock to ^20 shock terms of the (angle-dependent) 
shock radius, and the pre- and post-shock velocities [v-p, Vg) and densities (pp, ps), 



)E2 

^20, shock ' 



32^3/2 ^ 
715 dt, 

32^3/2 d 

15 At 



r-sh(e) 



j sin 9 j pv,- (3cos2 6l- l)dr+ J r^pVr{3cos^9~l)dr 



r^^{9) sin e'(3 cos^ 9-1) (ps^r.s - PpVr.p) <i9. 



d9 



(Bl) 



By exploiting the Rankine-Hugoniot jump conditions, the difference of the post-shock and pre-shock mass fiux can be 
eliminated and the shock velocity Vsh appears instead: 



Thus, we obtain 



\E2 

^20, shock 



327r3/2 d 



15 dt 



(Ps - Pp) Vsh = (Ps - Pp) hh- 

ipM-Ppmrli9)rli9){:icosH-l)sm9d9. 



(B2) 
(B3) 







'"sh(^) can be decomposed into spherical harmonics, or, in the case of axisymmetry, into Legendre polynomials, 

rsh{9) =Y.Pe{cos9)ae. (B4) 

Assuming a power-law dependence of the pre-shock density p on r'^ and a fixed ratio /3 of the post-shock and pre-shock 
density, we can express the angle-dependence of ps{9) and Pp{9) in terms of ae as well: 



Ps{9) - p,{9) ^{/3- l)p-p ( ^-^] ' ^{13- l)pp 



l-f7^P,( 



cosyja^ag -f 




)aeao ^ 



(B5) 

Here, p is an appropriate average value of p{9), and the ratio (?'sh(^)/ao)''' has been expanded into a power series. The 
resulting expression for ^fo shock '^^^^ expanded into a power series in the higher multipoles ai, 02, . . . of the shock 
position. Retaining only the linear terms and converting back to non-geometrized units, we finally obtain: 



2567^3/2 

^20,shock ~ . Pp ifi - 1) «o [(4 + 7) a2ao + a2ao] 



5V15 



(B6) 



^ For finite resolution, the identity J t^^ dV = J x^r^^dV = 
J r'^'x'x^ dV does not hold exactly due to discretization errors, 
but in practice the wave signal is hardly affected by this problem. 



Non-radial velocities occur behind a non-spherical shock, but 
a more detailed analysis of the jump conditions for oblique shocks 
reveals that this leads only to higher-order corrections to the re- 
sulting gravitational wave amplitude. 
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C. LEDOUX CRITERION, BRUNT-VAISALA AND PLUME FREQUENCY IN GENERAL 

RELATIVITY 

It may not be immediately obvious to the reader how buoyancy effects are to be treated in GR, and how, e.g., the 
analysis of the deceleration of convective downdrafts at the boundary of the convection zone presented bv lMurphv et all 
carries over to the relativistic case. We therefore briefly recapitulate the principles governing the motion of 
buoyancy-driven bubbles in the framework of GR hydrodynamics. Although the problem of convective stability has 
been studied in general relativity (Thornc 1966), we find it advisable to explicitly derive the equations for the particular 
gauge and frame-of-reference used in our paper in order to eliminate a possibl e source of conf usion. 

We start from the relativistic momentum equation in the 3 + 1 formulation (jBanvuls et al.l[l997i) . 

d^phW^, , dV^ [phW'v. (v^ - + SiP\ dg,^ 



dt dxi 2 dx3 

Here, p, h, W, P, Vi, T'"^ denote the rest-mass density, the specific enthalpy, the Lorentz factor, the pressure, the 
three- velocity and the stress-energy-tensor of the fluid, respectively. The four-metric is given in terms of the lapse 
function a, the shift vector /3% the conformal factor <p, and the flat-space three metric ^ij by Equation (jAip . which is 
repeated here for convenience: 



Finally, 7 and g denote the determinants of the three- metric 7^^ and the four- metric g/^,^^. 

We now consider the equation of motion for a convective blob of matter in a hydrostatic background medium which 
is displaced by a short distance Sr from its original position and maintains pressure equilibrium with its surroundings. 
The acceleration of the blob due to pressure gradients is given by the second term on the LHS of Equation (|Cip . 
which can be replaced by the gravitational source term for the background medium using the condition of hydrostatic 
equilibrium. The equation of motion (including the gravitational source term) is therefore given by 

where the subscripts denote whether T^'^ is evaluated for the blob or the background medium ("bg"). In order to 
simplify this result, we assume slow subsonic motion {v <^ c^), and retain density and energy perturbations only in 
the buoyancy term on the RHS (the Boussinesq approximation). The only remaining contribution on the RHS then 
contains T^jqi^ — T^g, and we obtain 

ph — ^-S{p + pe) — , (C4) 

where S{p + pe) denotes the (energy-)density contrast between the blob and the background medium (with e denoting 
the specific internal energy). For a small radial displacement Sr from the initial position, 5ip + pe) can be expressed 
as CL^r, Cl being the (relativistic) Ledoux criterion (cp. the "Schwarzschild discriminant" of lThornill966[ ). 

^^_ dpil + e) / dp(l + e) \ dP _ dpjl + 6) IdP 

^ dr \ dP J .^Y^^^^nst. 9r dr dr ' 

The equation of motion for the blob can thus be written as 

di^r,biob _ 0~^Cl da 
dt ph dr 

°^ .. ,_2dwr,biob aC-L da 

r = ad) = 7 — 

^ dt ph<j>'^ dr 

The blob therefore oscillates with an angular frequency N given by 

2 oCl da 



Sr, (C6) 
Sr. (C7) 



This is the relativistic Brunt- Vaisala-frequency. It should be noted that we have measured N in coordinate time and 
not in the frame of an observer comoving with the fluid. As the spacetime is stationary to very good approximation 
in our case, this implies that N is also the angular frequency seen by an observer at infinity^ 

For determining the "plume frequency" and the penetration depth Dp used in the model of iMurphv et al.l (|2009t ). 
we estimate the turnaround time T and Dp for a convective blob overshooting into a stable layer assuming a constant 
value of N equal to that at the turning point during the decelerating process. Moreover, we assume that the initial 
density contrast to the background medium at the boundary of the convective and non-convective regions is negligible 
compared to the density contrast after penetration into the non-convective region. Otherwise an additional acceleration 

* This can easily be seen by considering geodesies starting from , , r j ■ i u t-u u 1 ■ i- * 1 j-- jj- 

r ji I I I J / Trr j i- j , i oG transiormed into each other by applying a time translation 01 

the position of the blob at two different coordinate tmes t and everywhere. An observer at infinity therefore measures the same 

t + dt. Because of the stationarity of the metric, the geodesies can time difference 5t 
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term would have to be included in the equation of motion. Equation (jC7p can then be solved using 6r{0) — and 
r(0) — ^acjy^Vj. {vr being the plume velocity in the convective region) as initial conditions, 

Sr{t) = ac/j^Vr sin Nt. (C9) 

The penetration depth is therefore given by 

Dp = a(t)^vr, (CIO) 

and the characteristic (ordinary) frequency fp for the motion of the decelerating blob corresponding to the angular 
frequency is just at 

U = ^- (cii) 



